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We calculate the fluxes and energy spectra of high-energy (E > 50 GeV) neutrinos from the an- 
nihilations of supermassive (10 8 GeV < M < 10 16 GeV), strongly interacting dark matter particles 
in the core of the Sun. We take all significant aspects of neutrino propagation through matter into 
account, as well as oscillations in matter and vacuum. We also calculate the resulting event rates 
in an idealized 1 km 3 ice detector. We find that the signal should be well above background and 
easily observed by next-generation neutrino detectors such as IceCube. 
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PACS numbers: 13.15.+g, 14.60.Pq, 95.35.+d 



I. INTRODUCTION 

The identity of the non-baryonic dark matter, believed 
to account for about 30% of the total density of the Uni- 
verse jjj, is one of the major unanswered questions in 
cosmology today. There is much speculation that the 
dark matter consists of massive non-Standard Model el- 
ementary particles. The so-called "wimps," or weakly 
interacting massive particles, are the best known such 
candidate. These are often taken as the lightest super- 
symmetric particle, with a mass of no more than about 
7 TeV [^j . Any dark matter particle which is a thermal 
relic cannot have a mass of more than about 340 TeV || . 
It has been shown that wimps could be captured by the 
Sun [Q U , and that their subsequent annihilations could 
produce an observable flux of high-energy neutrinos || . 

Recently, another candidate for particle dark matter 
has been proposed, the "wimpzilla" [0, ||, [| |l^, [Tl| . 
These are gravitationally produced towards the end of 
inflation by the interaction of the wimpzilla field with the 
inflating space-time, and in sufficient abundance to be the 
dark matter. The wimpzilla is usually not assumed to be 
coupled to any other fields, although the case of coupling 
between the wimpzilla and the inflaton has been studied 
and wimpzilla production still found to be robust ; in- 
flaton decays have also been shown to be a possible source 
of supermassive particles [[l3| [lj] . The original wimpzilla 
calculations were done in the context of chaotic inflation 
(V(<p) = \m 2 ^(j) 2 ), but they have been shown to be abun- 
dantly produced for other inflaton potentials too Jl5| . 
Wimpzillas are most efficiently produced at extremely 
high masses, on the order of 10 12 GeV. 

A crucial difference between wimps and wimpzillas is 
that the latter are never in thermal equilibrium, and 
as such, their masses are not thermodynamically con- 
strained. It also follows from this that their present-day 
abundance does not depend on whether they have strong, 
weak, electromagnetic, or only gravitational interactions 
(although there are other considerations which generally 
rule out charged dark matter jl6|, |l7j|). 

We assume in this paper that the non-baryonic dark 
matter consists of wimpzillas with strong interactions, 
referred to from now on as "simpzillas." Simpzillas, like 



thermal wimps, can be captured by the Sun, and their 
annihilations in the solar core can produce high-energy 
neutrinos. We present the results of a detailed Monte 
Carlo calculation simulating the propagation of high- 
energy neutrinos from simpzilla annihilations through the 
Sun, and also taking oscillations in matter and vacuum 
into account. 

Our general conclusion is that the flux of high-energy 
neutrinos from simpzilla annihilations should be well 
above background for a broad range of parameter space, 
and observable by next-generation neutrino detectors 
such as IceCube. 



II. INITIAL NEUTRINO FLUX 

Albuquerque, Hui and Kolb Q, henceforth AHK, de- 
rived the capture rate of simpzillas by the Sun. The cap- 
ture of thermal wimps by the Sun had previously been 
studied in j|, ||. AHK assumed that the dark matter 
consists entirely of simpzillas in a Maxwell-Boltzmann 
velocity distribution, and that the interaction cross sec- 
tion a is of the order of the strong force, implying that 
the Sun is many interaction lengths thick. AHK's for- 
mulas for the capture rate are given in equations (2.5) 
and (2.7) of their paper. We plot the capture rate Tq 
in Figure |l|. The capture rate has two different forms, 
determined by the efficiency of the simpzilla energy loss 
in the Sun. 

Once captured by the Sun, the simpzillas rapidly fall 
to the core and annihilate with each other. AHK showed 
that equilibrium between annihilation and capture is 
reached very early in the lifetime of the Sun, and would 
obtain today. In equilibrium, by definition, Ta — Tc/2 
(each annihilation destroys two simpzillas). Note that 
this has the effect of making the equilibrium value of Ta 
dependent only on the interaction cross section a and not 
the annihilation cross section a a- 

High-energy neutrinos are produced by the simpzilla 
annihilations, which produce a quark or gluon pair that 
then fragment into hadronic jets containing a large num- 
ber of particles. AHK used the fragmentation function 
formalism of |lSfl to calculate the numbers of hadrons pro- 
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duced per annihilation. In the dense solar core, hadrons 
composed of light and charmed quarks lose most of their 
energy before decaying. Hadrons with bottom and top 
quarks, however, have much shorter lifetimes and de- 
cay before substantial energy losses. The neutrinos from 
these decays are at high energies. 

In this paper, as AHK did, we consider the high-energy 
neutrinos from top hadrons produced in simpzilla anni- 
hilations. Although high-energy neutrinos can also come 
from bottom hadrons, it is generally quite difficult to 
calculate their flux due to the large number of B me- 
son decay modes, the products of which may interact 
in the Sun before themselves decaying. Since we find 
that the neutrino flux from top hadron decays alone is 
well above background, including the neutrinos from B 
decays will only add to an already detectable signal, al- 
though most of these will be below the 50 GeV cutoff 
energy we choose. 

Approximately 2.8 x 10 5 V^i2 top hadrons are 
produced per simpzilla annihilation, where M\% = 
mx/W 12 GeV. Top quarks almost always decay in the 
channel t — > Wb. The W, in turn, decays with equal 
branching ratios (each about 10.5%) into ev e , and 
tv t . The W, like the t, decays before virtually any en- 
ergy loss. When the tv t pair is produced by the W decay, 
the r also decays before losing much energy, producing 
a second v T . About 18% of the time, it also produces a 
second v e and another 18% of the time a second v^. 

AHK showed that the top hadrons have an energy dis- 
tribution proportional to E^l 2 . They also calculated 
the distribution of the subsequent neutrinos. The total 
number of neutrinos produced above 50 GeV from this 
i — > W decay chain is 



dt 



(2.1) 



where £ denotes the neutrino flavor (e, [A, or r). m = 1 
for I = t and 1/2 for the other two flavors, representing 
the fact that roughly twice as many v T are produced per 
annihilation. 

The initial neutrino energy distribution calculated by 
AHK is 



= KdO 4 ^/M^T A ^E~ x 0.939 
dhjdt 

E + mw 



[E + mt] (E + mtf-m 2 (E + m w ) 



' w 



x Q(E-E min ) 



(2.2) 



50 GeV. This is near the lower limit of 
energies detectable by IceCube p0| , and also roughly 
the lowest neutrino energy possible in the two-body W 
decay. 

In Figure |^, we show the initial simpzilla neutrino flux 
at the core of the Sun. 
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FIG. 1: The capture rate of simpzillas by the Sun as a 
function of simpzilla mass mx, for three different choices of 
the interaction cross section a. In the shallower parts of the 
curves, most of the simpzillas are captured. In the steeper 
parts, the simpzilla energy loss in the Sun is inefficient and 
only those with lower velocities are captured. 
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FIG. 2: The initial energy distribution of neutrinos from 
top hadrons produced in simpzilla annihilations. Kg — 1 for 
v T and 1/2 for and v^. Note that the typical neutrino 
energy is much lower than the simpzilla mass, due to the 
large numbers of neutrinos produced per annihilation. 



III. NEUTRINO INTERACTIONS 



Neutrinos with energies above about 100 GeV have 
significant interactions with matter as they propagate 
through the Sun. The majority are charged-current (CC) 
interactions with nucleons. The rest are neutral-current 
(NC) scatterings; the cross sections for the latter are 
about one-third of the charged-current scattering cross 
sections at the same energy. Neutrino-electron scat- 
terings in this energy range are negligible compared to 
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yl , and accordingly we do 



neutrino-nucleon scatterings 1 
not consider them in our code. 

We have calculated the charged and neutral-current 
neutrino-nucleon cross sections between 1 GeV and 
10 9 GeV. The most recent previous calculations in this 
energy range 21, have assumed an isoscalar nucleon; 
that is, one in which the quark distribution functions are 
the average of those for the proton and neutron. This 
is a good approximation for a medium such as rock (the 
calculations were in the context of high-energy neutrino 
beams propagating through the Earth). However, the 
proton number density in the Sun is between two and 
six times the neutron number density, and the isoscalar 
approximation does not hold. We show the ratio of pro- 
ton and neutron number densities in Figure [j[ calculated 



using the Standard Solar Model of 23 



We have used the recently published CTEQ6-L (lead- 
ing order) parton distribution functions |24j, together 
with the deep inelastic scattering formalism of fl25fl , to 
calculate the cross sections. Given the energies involved, 
we may neglect the electron and muon masses, and so the 
v e and CC cross sections are identical. The r mass, 
however, significantly suppresses the v T CC cross sec- 
tions up to about 50 GeV. Although the effect of this on 
our calculation of the emergent simpzilla neutrino fluxes 
is minor (we only consider energies above 50 GeV, and 
neutrino interactions in the Sun are significant only above 
about 100 GeV), we have taken it into account. The NC 
cross sections are identical for all three flavors. We show 
our results in Figures 

Note that for E v < 10 4 GeV (and in the case of u T , 
above energies where the r mass suppression is signifi- 
cant), the cross sections are approximately proportional 
to energy. At higher energies, the rate at which the cross 
section increases is suppressed by the gauge boson prop- 
agator pi]] . However, almost all of the simpzilla neutrino 
flux is below 10 4 GeV. 

We have also calculated the differential cross sections 
da/dy, where the inelasticity parameter y is 



E' 

V = i 



(3.1) 



and E' is the final charged lepton (neutrino) energy in 
a CC (NC) scattering. We show these for v e /v^ and v T 
CC scattering in Figures ^ and ||; the others are similar. 
Our code uses these differential cross sections to obtain 
probability distributions for E' . 

We find numerically that almost all neutrino interac- 
tions take place in the deepest part of the Sun, where the 
density is greatest, out to about O.li? . This is signifi- 
cant when oscillations are taken into account, because at 
high energies the oscillations and interactions effectively 
decouple in radius (see section V). 




FIG. 3: The ratio of the proton and neutron number densities 
in the Sun as a function of radius. The number of neutrons is 
greatest in the core where more hydrogen has been processed 
into helium. Isotopes other than 1 H and 4 He are neglected. 
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1 Except for v £ in the vicinity of the Glashow resonance at 6.3 PeV. 
However, the fluxes we consider are far below this. 



FIG. 4: The charged-current interaction cross section be- 
tween an electron or muon neutrino and a neutron (n), proton 
(p), and isoscalar nucleon (N). 



IV. CHARGED LEPTON INTERACTIONS 

The charged leptons produced in charged-current in- 
teractions are very important in determining the emer- 
gent neutrino flux. Electrons and muons tend to be 
stopped by their electromagnetic energy losses. How- 
ever, taus have very short lifetimes. They tend to decay 
before losing any significant fraction of their energy, and 
so the resulting v r are also high-energy. CC interactions, 
in effect, absorb the v e and but regenerate the v T . 
This phenomenon was first remarked upon by Ritz and 
Seckel || in the context of a calculation of neutrino spec- 
tra from thermal wimp annihilations in the Sun (see also 
p6|); Halzen and Saltzberg |27j have shown it to obtain 
for high-energy neutrino beams propagating through the 
Earth (see also ||, ||, |(J). 

In the absence of oscillations, this means that high- 
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FIG. 5: The charged-current cross sections for a tau neutrino. 
The electron/muon neutrino cross sections are also shown to 
illustrate the kinematic suppression due to the r mass. 
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FIG. 6: The neutral-current cross sections for all three fla- 
vors. 



energy v r are much more likely than high-energy v e and 
to survive. The v T are moderated down to lower ener- 
gies, to the point where the probability of further interac- 
tions is small. However, the total number of v T emerging 
is essentially equal to the initial number, whereas the v e 
and are attenuated. We will discuss the properties of 
the v T versus v^jv^ emergent fluxes in greater detail in 
the penultimate section, including the effects of oscilla- 
tions. 

In this section, we consider the energy losses and de- 
cays of the charged leptons in the solar medium. 



A. Energy losses 

As the charged leptons move through the Sun, they 
have electromagnetic interactions with the medium and 
lose energy. (Weak interactions, apart from decays, 
are insignificant for charged leptons at energies below 
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FIG. 7: CC differential cross sections for a v e or f M . From 
top to bottom at y = 5 x 10 -3 , the curves correspond to 
10 GeV, 10 2 GeV, 10 9 GeV. The cutoffs in the 10 GeV 
curves at y m 3 x 10 -3 are where the Q 2 of the gauge boson 
decreases below Aq CD , and hence the parton formalism used 
to calculate deep inelastic scattering becomes invalid. 
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FIG. 8: CC differential cross sections for a u T . The energies 
are the same as in Figure 0. The cutoffs at large y in the 
10 GeV and 100 GeV curves are due to the r mass. 

10 16 GeV Iplj ). The general equation for energy loss as 
a function of column depth x is 

-(f)- + E^, («) 

where a represents interactions with electrons and {/?■'} 
represent interactions with nuclei, with the summa- 
tion being over the different processes: bremsstrahlung, 
pair production, and photonuclear scattering. Below 
10 4 GeV, where most of the simpzilla neutrinos (and 
charged leptons from their CC scatterings) are produced, 
electronic losses dominate. In most terrestrial situations, 
these are ionization losses which can be calculated with 
the Bethe-Bloch formula |3^, |33| , but inside the Sun the 
electrons are unbound. Therefore, a is the energy loss 
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rate of a charged particle in a plasma 



a = (9.2 x 10 6 



Z 



Al3 2 c 



In (2m e 7/3 2 ) - In 



(4.2) 



The units of a in Eq. fly) are GeV/(g/cm 2 ). The 
quantities Z and A are the atomic number and weight 
of the nuclei in the medium - for the Sun, which is (ap- 
proximately) a mixture of 1 H and 4 He, we calculate the 
energy losses separately for the two isotopes and then 
weight them by their mass fractions. The (3 and 7 are 
the Lorentz parameters of the charged lepton and n e is 
the electron number density. 

The nuclear energy losses become dominant at energies 
above about 10 4 GeV, although as mentioned they are 
not very important for the simpzilla neutrino flux which 
tends to be at energies below this. The formulas for {Z?- 7 } 
are given, for example, in pl| . 

Our code simulates both electronic and nuclear energy 
losses for /x and r. We do not follow e since they never 
decay. In Figure ^, we show the muon and tau ranges and 
decay lengths in the core of the Sun, which is where the 
majority of neutrino interactions take place. The range 
is defined as the distance, in the absence of decays, over 
which the lepton would lose all its energy; and similarly 
the decay length is the (mean) distance, in the absence 
of energy losses, the lepton would travel before decaying. 



B. Decays 

Both muons and taus decay back into neutrinos, but 
the muons have lost almost all their energy by this point 
and essentially decay at rest. The resulting are below 
50 MeV and do not concern us. The taus, on the other 
hand, have negligible energy losses except at very high 
energies (E > 10 7 GeV), which again does not concern 
us since very little of the initial simpzilla neutrino flux is 
in this range. The v T from the decay carries on average 
about 2/5 the energy of the r, or about 1/4 the energy 
of the incident v r (the mean energy of the r relative 
to the incident v T can be obtained from the inelasticity 
parameter distributions; in the energy range we consider, 
the mean value of y is about 0.4, higher than for the 
energies at which neutrinos interact significantly in the 
Earth). 

One phenomenon our code takes into account is the 
production of secondary i> e and in r decays, as dis- 
cussed by J3(| . Approximately 18% of r decays produce 
a (for r^, a (vy)), and another 18% produce a v e . 
These secondary neutrinos are created with about half 
the energy of the accompanying z/ T , or about 1/8 the 
energy of the incident v T . When the v T experiences its 
last few charged-current scatterings and is downscattered 
to energies at which it escapes the Sun, these secondary 
neutrinos emerge too. Ref. |3C|] considered high-energy 
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FIG. 9: Decay lengths and ranges of fi and r in the solar core 
(where most charged-current interactions take place). Note 
that except at the highest energies, the r decay length is much 
shorter than its range; hence it decays before losing virtually 
any energy. The opposite is true for the /x, which is essentially 
brought to rest before decaying. 



neutrino beams from extragalactic point sources prop- 
agating through the Earth, and showed that the addi- 
tional flux of secondaries could substantially modify the 
detected signal. 

The major decay channels of r are r — > v T [LVy (18%); 
r -» v T eu e (18%); r -> v t t: (12%); r -> v T p (26%); 
t — > v T a\ (13%); and r — > v T X (13%), where X indicates 
other massive hadrons. The combined v T decay distribu- 
tion is given in ||l| . We show it in Figure [l0[ along with 
that of the secondary v e and which is straightforward 
to obtain (see e.g. Q). Our code randomly samples over 
the appropriate distributions to obtain the energies of the 
neutrinos from r decays. 



V. NEUTRINO OSCILLATIONS 

Both matter and vacuum oscillations are significant for 
simpzilla neutrinos, and have a large effect on the emer- 
gent flux. Our code takes the oscillations of all three 
flavors into account. We use values for the mass-squared 
differences and mixing angles which are consistent with 
recent observations. From the SuperKamiokande results 
p| [3§, we choose Am^ = 3 x 10~ 3 eV 2 , sin 2 £ = 0.1 
(which is the upper limit on £), and sin 2 9 = 0.5, where 
the angles in the three-flavor mixing matrix are as in Eq. 
(3.1) of |37]]. The other parameters are constrained by so- 
lar neutrino observations, such as those recently recorded 
by the Sudbury Neutrino Observatory (SNO) g§|, @. 
Analysis of the SNO data favors the "large mixing angle" 
(LMA) solution pd), and so we choose values character- 
istic of this: Am^7 = 2 x 10~ 5 eV 2 and sin 2 lo = 0.2. We 
assume the normal mass hierarchy with m\ < mi < m^, 
and set the CP-violating phase in the mixing matrix to 
0. 
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FIG. 10: Distributions of neutrino energies from r decays, 
with z = Ev/Et and the r assumed to be ultrarelativistic. 
Solid curve: v r energy distribution, summing over all signifi- 
cant decay channels of the r and weighting by the branching 
ratios. Some of the channels produce massive hadrons and 
become kinematically impossible above a certain z, the rea- 
son for the jagged appearance. Dashed curve: secondary 
and v e decay distributions from the 18% of r decays which 
produce them. Dotted curve: the v T energy distribution in 
the secondary-producing channels. We assume the r to have 
the parity it would have if produced in a v T CC interaction 
(for t~, negative; for r + , positive). 

To simulate the oscillations, we evolve the flavor ampli- 
tudes using the analytical solution for the time-evolution 
operator in jl2|, 0. Our Monte Carlo step sizes are 
sufficiently small that the oscillation probabilities do not 
change significantly over the step, and the density at each 
step may be taken as constant. If the neutrino is deter- 
mined to have had a charged-current interaction during 
the step, we sample over the current probabilities to find 
the flavor of the resulting charged lepton; if the lepton 
subsequently decays at high energy, the resulting neu- 
trino begins purely in that flavor. 

In Figures [ll] and |lj, we show the oscillation proba- 
bilities of 100 GeV neutrinos and antineutrinos traveling 
through the Sun without interacting (so their energies 
do not change). This is roughly the lower limit of where 
neutrino interactions in the Sun are significant. 

We have noted above that most neutrino interactions 
take place between the center of the Sun and O.li?©. At 
higher energies, the oscillation wavelengths increase, to 
the point where the oscillation probabilities do not vary 
substantially over this distance. Interactions at energies 
much greater than 100 GeV thus effectively "decouple" 
in radius from oscillations (though if our calculation were 
extended to include sterile neutrinos, this might not hold; 
see @). 

We also would not expect different reasonable choices 
of the oscillation parameters to make a difference of more 
than a factor of two or so in the detection rates. Real- 
world detectors most easily observe the muons created 
in and charged-current interactions (the v e and v T 



FIG. 11: Oscillation of 100 GeV neutrinos propagating 
(without interactions) from the center of the Sun to its sur- 
face. From top to bottom: oscillation probabilities of 
and v T , Note that the oscillations of u e are almost completely 
suppressed until the resonance. 

signals are much less directional) , so the relevant issue is 
how much different parameter choices would change the 
emergent + flux. As can be seen from Figures [ll] 
and [121 the worst-case scenario would completely deplete 
the flux by maximizing the v e conversion probabilities; 
however, the flux would be unaffected since the v e do 
not encounter a resonance and barely oscillate. (A dif- 
ference in the sign of Ato^, which is currently unknown 
and which we assume to be positive, would cause the v e 
rather than the v e to have a resonance, so the graphs 
would essentially be interchanged.) 

The <-» v T oscillations in the Sun are vacuum oscil- 
lations with the atmospheric mixing parameters. The v e 
resonance is governed by the mixing angle £ (correspond- 
ing to Uei) and Am^. This is in contrast to normal solar 
neutrinos, whose MSW resonance is governed by Am^ 
and u>. For the simpzilla neutrinos, that resonance is in 
the outermost layers of the Sun at very low densities, 
and is not numerically significant compared to the other 
resonance. We note that vacuum oscillations outside the 
Sun are completely averaged out by the detector energy 
resolution (see section VI). We also note that although 
the detector is considered to be looking at upcoming sim- 
pzilla neutrinos which have come through the Earth, the 
density of the Earth is too high for MSW oscillations 
inside it to be important. 



VI. NEUTRINO FLUX AT EARTH 

U sing the initial neutrino energy distribution given in 
( J2.2D and taking all the physical effects discussed above 
into account, we have calculated the resulting neutrino 
flux at Earth. We have also, with IceCube in mind, cal- 
culated the event rates it would produce in an idealized 
1 km 3 ice detector. 



7 



100 GeV gg, F K , V T 




FIG. 12: Oscillations of 100 GeV antineutrinos propagating 
from the center of the Sun to its surface. The v R , unlike the 
v e , does not encounter a resonance, and its oscillations are 
minimal. 



A. Monte Carlo code 

We have written a Monte Carlo code to simulate the 
propagation of high energy neutrinos through matter. 
Each neutrino is followed as it passes through the Sun. 
The step sizes are chosen to be much smaller than both 
the interaction length and the neutrino oscillation length. 

At each step, we use a pseudo-random number gener- 
ator to determine whether an interaction occurred over 
that step, whether it was a charged- or neutral-current 
interaction, and whether it occurred on a proton or neu- 
tron. In the event of a neutral-current interaction, the fi- 
nal neutrino energy is determined by randomly sampling 
over the appropriate inelasticity parameter distribution 
(see section III). 

Charged-current interactions are effectively measure- 
ments of the neutrino flavor. As such, when a charged- 
current interaction occurs, we randomly sample over the 
current oscillation probabilities to determine the flavor 
of the neutrino and outgoing lepton. When and if the 
lepton decays back into a neutrino, it starts in a pure 
flavor state. We calculate the initial energy of the lepton 
using the inelasticity parameter distribution. 

The fi and r produced in charged-current interactions 
are followed until they decay or their energies fall below 
50 GeV. The e produced in v e charged-current interac- 
tions are not followed since they never decay. (In prac- 
tice, the n almost never decay at high energies either, 
but we follow them for consistency.) We calc ulate their 
energy losses over each step according to (4.1). We also 
use their boosted decay lengths to randomly determine 
whether a decay occurs on each step. If so, we sample 
over the appropriate decay distribution to determine the 
energy of the outgoing neutrino. We also calculate the 



energies of the secondary v e and in r decays using the 
energy distributions for those modes. 

Neutrino oscillations are simulated both in the Sun 



and in the vacuum between the Sun and the Earth. We 
average the oscillation probabilities over 10% of the Sun- 
Earth distance, which is mathematically equivalent to a 
10% energy resolution in the detector (characteristic of 
real- world detectors) . 

The detector is considered to be able to detect every 
neutrino charged-current interaction above 50 GeV tak- 
ing place in a cubic kilometer of ice. In practice, real 
detectors are considerably more complex than this and a 
sophisticated treatment would require a detector Monte 
Carlo, which is beyond the scope of this work. The de- 
tector is also assumed to have an angular resolution of 
1°, which is approximately the size of the Sun on the sky. 

Having calculated the neutrino flux at the Earth, we 
calculate the detection rate by integrating the product 
of the flux, the area of the detector (1 km 2 ), and the 
probability of the neutrino having a charged-current in- 
teraction, the integration being over the energy E. The 
probability of a charged-current interaction in the detec- 
tor is 



P cc (E) = l-e 
with the interaction length 

Lint{E) 



-X/L int (E) 



(6.1) 



(6.2) 



and A = 1 km is the size of the detector. The number 
densities n n and n p are for ice, and the charged-current 
cross sections a cc (E) were discussed in section III. 

The background flux from atmospheric neutrinos 
above 50 GeV is approximately J45[ 



d<$>atr, 



dE dA dt duj 
[km" 2 



= (1.1 x 10 12 )£~ 3 - 2 
1 deg" 2 GeV" 1 ] . 



(6.3) 



Integrating this over energy, the detection probability, 
the detector area, and the angular size of the Sun gives 
about 2 atmospheric neutrino events per year above 
50 GeV which come from the direction of the Sun. These 
will be almost exclusively muon neutrinos, if we assume 
the detector is looking at upcoming neutrinos. We take 
this as our background. 

We also note that the Baksan neutrino telescope has 
placed a 90% confidence level limit of muon fluxes from 
nonatmospheric neutrinos coming from the direction of 
the Sun of about 10 4 km -2 yr _1 We take this as a 
rough upper limit. The future analysis of the data from 
AMANDA |47) may lower this limit. 



B. Results 



In Figures 13 and [16|, we show our numerical results 
for the simpzilla neutrino flux at the Earth. In Figures 
|i"7| and [l8l we show the event rates as a function of sim- 
pzilla mass mx for three different choices of the interac- 
tion cross section a. We show two cases: when neither 
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oscillations nor secondaries are included, and when both 
are included. We have found that the secondaries alone 
do not make a dramatic difference to the flux, largely be- 
cause it falls over several decades in energy, and much of 
it is in a region where neutrino interactions are not sig- 
nificant and thus secondaries are not produced. We ac- 
cordingly do not show the other two permutations. Pre- 
viously, only the case without oscillations or secondaries 
has been studied @ ff|] . 

We discuss this case briefly. As can be seen from Figure 
|l3| , the emergent v T fluxes exceed those of the other two 
flavors by more than the initial factor of 2 (see section 
II). This is due to the fact that the v T which experience 
charged-current interactions are not, as are the v e and 
v^, absorbed, but rather regenerated at lower energies, 
as remarked upon in section IV. 



The unscattered components of the v e 
beams relate to their initial fluxes as fllql 



and 



-E/E k 



dEdt dEdt 



(6.4) 



where is the "transparency energy," defined as the 
energy such that the mean number of interactions ex- 
perienced by the neutrino is 1. Equation (5.4) obtains 
because the cross sections are approximately linear func- 
tions of energy in the range we are considering. Numeri- 
cally, we find that E). ~ 130 GeV for v e and v^, 160 GeV 
for v Tl 200 GeV for D e and P„, and 230 GeV for v T . We 
use ( |6.4[ ) as the emergent flux for the electron and muon 
neutrinos 2 . 

To model the scattered component of the v r , we use the 
observation of ourselves and several other groups that the 
scattered v r emerge in a roughly log-normal distribution. 
We find numerically that about 80% of the emergent v T 
and v T are scattered. The log-normal distribution is de- 
fined as 



dn 
dE 



1 



^InlOaE 



exp 



1 

"2a 



2 1^ 



(6.5) 



and when transformed to a distribution in log-E, is a 
Gaussian with mean log E t and standard deviation a 
(whose units are decades). Note that E t ^ Ek, although 
they turn out to be similar. We show how the numeri- 
cal results compare to a log-normal distribution in Fig- 
ure [l4|. The correspondence is not exact, and the log- 
normal fit overestimates the high-energy neutrinos while 
underestimating the low-energy ones. However, it is ev- 
idently a decent approximation. The v T log-normal fit 
has a = 0.53 and E t — 60 GeV; the v T has a — 0.49 and 
Et = 113 GeV. 



2 A small number of u e and (as well as u T ) experience repeated 
neutral-current scatterings without any charged-current scatter- 
ings. We neglect this neutral-current-scattered component here, 
although it is taken into account by our code. 



Without secondaries or oscillations 




Energy [GeV] 

FIG. 13: The simpzilla neutrino fluxes at the Earth, neglect- 
ing both neutrino oscillations and secondary neutrinos from 
r decays. 



In Figure [15], we plot these analytical approximations. 
The v e and fluxes are in qualitative agreement with 
the numerical results. The v T fluxes are somewhat lower 
at low energies and higher at high energies than are the 
numerical results, again due to the log-normal underes- 
timating the fi rst while overestimating the second. We 
note that Eq. differs from Eq. (4.6) of |J], 

leading 

to a much smaller difference between the v T and v^jv^ 
event rates; this is due to our correctly taking the 1/E 
term in the Jacobian factor in the log-normal distribution 
into account. 

The oscillations clearly have a substantial effect on 
the flux. The v r are depleted, whereas the Vn and v e 
are enhanced. In a real world detector, this is advanta- 
geous, because it is generally easiest to observe muons 
from charged-current interactions. The detectability 
of the simpzilla neutrino flux is therefore improved by 
oscillations. 

In Figure [l^, we show the excludable region of sim- 
pzilla parameter space; that is, the range of rrix and 
a which would be ruled out by the failure to observe a 
high-energy neutrino signal coming from the direction of 
the Sun. Even better, of course, would be the observa- 
tion of such a signal. However, a more detailed analy- 
sis, which we leave to others, would be necessary to rule 
out other possible sources such as thermal wimps, par- 
ticularly massive ones that could produce neutrinos in 
roughly the same energy range. 



VII. CONCLUSIONS 

We have shown that the annihilations of simpzillas in 
the Sun should produce a conspicuous high-energy neu- 
trino signal in a 1 km 3 ice detector for a fairly large range 
of simpzilla mass and interaction cross section. Much of 
this range of parameter space has not yet been ruled out 
by current observations. 
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v T , E = 10 3 GeV 

1.0 1 1 




Energy [GeV] 

FIG. 14: Distribution of emergent v T energies for a monoen- 
ergetic beam with initial energy E a = 10 3 GeV. We have 
calculated the emergent distribution for other choices of Eq 
too, and found it to be generally independent of Eq. We plot 
as a distribution in log E rather than E, and neglect oscilla- 
tions. The distribution corresponds approximately (although 
obviously not exactly) to the log-normal distribution with the 
parameters shown. 



W/o secondaries, oscillations (analytical) 




1 2 1 3 

Energy [GeV] 



FIG. 15: Analytical approximation for the emergent sim- 
pzilla neutrino fluxes, neglecting secondaries and oscillations. 
We assume the z^ e and fluxes to consist entirely of the 
remaining unscattered neutrinos, and the v T flux to be the 
sum of an unscattered component and a scattered component 
obeying a log-normal distribution. We make similar assump- 
tions about the antineutrinos. The parameters used are given 
in the text. 



Our flux calculation takes all important physical ef- 
fects into account, including neutrino oscillations, which 
approximately equalize the fluxes of the different flavors. 
In calculating the event rates, we have made some simpli- 
fying assumptions about the detector, and a more real- 
istic calculation (say, for IceCube) would take the actual 
detector response characteristics into account. However, 
we would not expect this to alter our general conclusion 
that the simpzilla neutrino signal would be observable in 



With secondaries and oscillations 




1 2 1 3 

Energy [GeV] 



FIG. 16: The simpzilla neutrino fluxes at the Earth cal- 
culated with our Monte Carlo, taking both oscillations and 
secondaries into account. 



Without secondaries or oscillations 
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FIG. 17: The simpzilla neutrino event rates in a 1 km 3 ide- 
alized detector, neglecting secondaries and oscillations. The 
rate is plotted as a function of the simpzilla mass mx- We 
show three different values for a, the simpzilla interaction 
cross section: from left to right, these are a = 1CP 26 cm 2 , 
1CF 24 cm 2 , and 1CP 22 cm 2 . The background from atmo- 
spheric neutrinos (y^ and v^) is approximately 2 yr -1 . An 
upper limit from Baksan is about 10 4 yr -1 . 



such a detector. 

The next generation of neutrino detectors should be 
able to either rule out a large range of simpzilla param- 
eter space, or provide observational evidence for them. 
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FIG. 18: The simpzilla neutrino event rates, taking both 
secondaries and oscillations into account. The same choices 
for a are plotted as in Figure [l?]. 




FIG. 19: The regions of simpzilla parameter space which 
produce neutrino fluxes (in each flavor) that are above and 
below the background event rate of 2 yr" 1 . If no high-energy 
neutrino signal is observed from the Sun, it would definitively 
exclude the range of simpzilla mass and cross section to the 
left of the line. The plot was generated using the event rates 
which include both secondaries and oscillations. 
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